ICRR-Report-315-94-10, TU-457 



Gravitino Production in the Inflationary Universe and the 
Effects on Big-Bang Nucleosynthesis 

M. Kawasaki 

Institute for Cosmic Ray Research, The University of Tokyo, Tanashi 188, Japan 

T. Moroi * 

Department of Physics, Tohoku University, Sendai 980, Japan 

(March 29, 1994) 



^J- Abstract 

Gravitino production and decay in the inflationary universe are reexanimed. 
Assuming that gravitino mainly decays into photon and photino, we have cal- 
^ culated the upperbound of the reheating temperature. Compared to previous 

works, we have essentially improved the following two points; (i) the helicity 
CN ±| gravitino production cross sections are calculated by using the full relevant 

terms in the supergravity lagrangian, and (ii) high energy photon spectrum is 
obtained by solving the Boltzmann equations numerically. Photo-dissociation 
of the light elements (D, T, 3 He, 4 He) leads to the most stringent upperbound 
£^ of the reheating temperature, which is given by (10 6 -10 9 )GeV for the grav- 

itino mass lOOGeV-lTeV. On the other hand, requiring that the present mass 
density of photino should be smaller than the critical density, we find that 
the reheating temperature have to be smaller than (10 11 -10 12 )GeV for the 
photino mass (10-100)GeV, irrespectively of the gravitino mass. The effect 
of other decay channel is also considered. 
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I. INTRODUCTION 



When one thinks of new physics beyond the standard model, supersymmetry (SUSY) 
is one of the most attractive candidates. Cancellation of quadratic divergences in SUSY 
models naturally explains the stability of the electroweak scale against radiative corrections 
[1,2]. Furthermore, if we assume the particle contents of the minimal SUSY standard model 
(MSSM), the three gauge coupling constants in the standard model meet at ~ 2 x 10 16 GeV 
[3,4], which strongly supports the grand unified theory (GUT). 

In spite of these strong motivations, no direct evidence for SUSY (especially superpart- 
ners) has been discovered yet. This means that the SUSY is broken in nature, if it exists. 
Although many efforts have been made to understand the origin of the SUSY breaking, we 
have not understood it yet. Nowadays, many people expect the existence of local SUSY (i.e., 
supergravity) and try to break it spontaneously in this framework. In the broken phase of 
the supergravity, super-Higgs effect occurs and gravitino, which is the superpartner of gravi- 
ton, acquires mass by absorbing the Nambu-Goldstone fermion associated with the SUSY 
breaking sector. In many broken local SUSY models, we expect that the mass of gravitino 
m3/2 is at the same order of those of squarks and sleptons since the following (tree level) 
super-trace formula among the mass matrices Ai 2 /s holds [5]; 

StiM 2 = ]T (-1) 2J (2 J + l)tiM 2 j ~ 2(n - l)m 2 /2 , (1) 

spin J 

where n is the number of chiral multiplet. For example in models with the minimal kinetic 
term, this is the case and all the SUSY breaking masses of squarks and sleptons are equal 
to gravitino mass at the Planck scale. But contrary to our theoretical interests, we have no 
hope to see gravitinos directly in the collider experiments since the interaction of gravitino 
is extremely weak. 

On the other hand, if we assume the standard big-bang cosmology, the mass of gravitino is 
severely constrained. If gravitino is unstable, it may decay after the big-bang nucleosynthesis 
(BBN) and produces an unacceptable amount of entropy, which conflicts with the predictions 
of BBN. In order to keep success of BBN, the gravitino mass should be larger than ~ lOTeV 
as Weinberg first pointed out [6]. Meanwhile, in the case of stable gravitino, its mass should 
be smaller than ~ IkeV not to overdose the universe [7]. Therefore, the gravitino mass 
between ~ IkeV and ~ lOTeV conflicts with the standard big-bang cosmology. 

However, if the universe went through inflation, we may avoid the above constraints 
[8] since the initial abundance of gravitino is diluted by the exponential expansion of the 
universe. But even if the initial gravitinos are diluted, the above problems still potentially 
exist since gravitinos are reproduced by scattering processes off the thermal radiation after 
the universe has been reheated [9-17]. The number density of secondary gravitino is propor- 
tional to the reheating temperature and hence, upperbound on the reheating temperature 
should be imposed not to overproduce gravitinos. Therefore, even assuming the inflation, 
a detailed analysis must be made to obtain the upperbound on the reheating temperature. 
The case of stable gravitino is analyzed in Refs. [13,15,17] and we will not deal with it. 

In this paper, we assume that the gravitinos are unstable and derive the upperbound 
on the reheating temperature. The analysis of the cosmological regeneration and decay of 
unstable gravitino has been done in many articles. These previous works show that the most 
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stringent upperbound on the reheating temperature comes from the photo-dissociation of 
the light nuclei (D, T, 3 He, 4 He). Once gravitinos are produced in the early universe, most 
of them decay after BBN since the lifetime of gravitino with mass O(100GeV — lOTeV) 
is O((10 8 — 10 2 )sec). If gravitinos decay radiatively, emitted high energy photons induce 
cascade processes and affect the result of BBN. Not to change the abundance of light nuclei, 
we must constrain the number density of gravitinos, and this constraint is translated into 
the upperbound of reheating temperature. 

In order to analyze the photo-dissociation process, we must calculate the following two 
quantities precisely; the number density of gravitinos produced after the reheated universe, 
and the high energy photon spectrum induced by radiative decay of gravitinos. But the 
previous estimations of these values are incomplete. As for the number density of gravitino, 
most of the previous works follow the result of Ref. [11], where the number density is 
underestimated by factor ~4. Furthermore, in many articles, the spectrum of high energy 
photon, which determines the photo-dissociation rates of light elements, are calculated by 
using a simple fitting formula. In this paper, we have treated these effects precisely and found 
the more stringent upperbound on the reheating temperature than the previous calculations. 
The plan of this paper is as follows. In Sec. II, we calculate the gravitino production cross 
section in the early universe. In Sec. Ill, high energy photon spectrum induced by the 
radiative decay of gravitino is obtained by solving the Boltzmann equation numerically. 
The results are shown in Sec.V. Other cosmological constraints is considered in Sec. VI and 
Sec. VII is devoted to discussions. 



II. GRAVITINO PRODUCTION IN THE EARLY UNIVERSE 

After the universe has reheated, gravitinos are reproduced by the scattering processes 
of the thermal radiations and decay with decay rate of order of m^/M^ where M p i — 
\/8tvM ~ 1.22 X 10 19 GeV is the Planck mass. Since the interaction of gravitino is very 
weak, gravitino cannot be thermalized if the reheating temperature is less than 0{M p i). 
In this case, Boltzmann equation for the gravitino number density n 3 /2 can be written as 

dn 3 / 2 , . 2 m 3/2 n 3/2 

—— + 3Hn 3/2 = {^totVrei} n rad - y— -r , (2) 

at (E 3/2 ) r 3/2 



where H is the Hubble constant, E ioi is the total cross section which is derived in appendix A, 
(• • •} means thermal average, n ra( i = ((3)T ?j /tt 2 represents the number density of the scalar 
boson in thermal bath, v re i is the relative velocity of the scattering radiations ((v re i) — 1 in 
our case), and m^^l ^3/2) is the averaged Lorenz factor. Note that the first term of right 
hand side (r.h.s.) of Eq.(2) represents contribution from the gravitino production process, 
and the second one comes from the decay of gravitino. In Eq.(2), we have omitted the 
terms which represents the inverse processes since their contributions are unimportant at 
low temperature that we are interested in. For the radiation dominated universe, H is given 
by 

D / AT 

(3) 
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where R is the scale factor and N* is the total number of effectively massless degrees of 
freedom. For the particle content of MSSM, N*(T R ) ~ 228.75 if T R is much larger than the 
masses of the superpartners, and N*(T <C lMeV) ~ 3.36. 

At the time right after the end of the reheating of the universe, the first term dominates 
the r.h.s. of Eq.(2) since gravitinos have been diluted by the de Sitter expansion of the 
universe. Using yield variable Y3/2 = ^3/2/^™^ an d ignoring the decay contributions, Eq.(2) 
becomes 

dT = HT ' 1 ' 

where we have assumed the relation 

RT — const. (5) 

Ignoring the small T-dependence of S ioi , we can solve Eq.(4) analytically. Integrating Eq.(4) 
from the reheating temperature T R to T (T R ^> T) and multiplying the dilution factor 
Ns(T)/Ns(Tr), the yield of gravitino is found to be 

v rr\ N s( T ) v n rad (T R ) {Z tot v ret } 

Y3,2lT) = w x mfl) ■ (6) 

For the MSSM particle content, N S (T R ) ~ 228.75 and N S (T < lMeV) ~ 3.91. Eq.(6) 
shows that Y3/2 is proportional to T R . Note that the numerical value of Y3/2 in our case is 
about 4-5 times larger than the result in Ref. [11]. Some comments on this difference are 
given in section VII. From Eq.(6), we can derive the simple fitting formula for Y3/2; 

>Wr<lMeV) = 2.14X10-" (^) {l - 0.0232 log (^) } , (7) 

where the logarithmic correction term comes from renormalization group flow of the gauge 
coupling constants. The difference between the exact formula (6) and the above approxi- 
mated one is within ~ 5% (~ 25%) for 10 6 GeV < T R < 10 14 GeV (10 2 GeV < T R < 10 19 
GeV). 

As the temperature of the universe drops and H^ 1 approaches r 3 / 2 , decay term becomes 
the dominant part of the r.h.s. of Eq.(2). Ignoring the scattering term, Eq.(2) can be 
rewritten as 



dt r 3/2 ' 



(8) 



where we have taken m 3 / 2 / ^-£-3/2) — 1 since gravitinos are almost at rest. Using Eq.(6) as 
a boundary condition, we can solve Eq.(8) and the answer is 

V (t\ - - V n ™d{T R ) {^totVrel) ( t \ 

n Ta d\t) N S {T R ) H(T R ) \ t 3/2 J 

where the relation between t and T can be obtained by solving Eq.(3) with Eq.(5); 



1 90M 2 „ . , 
t = -J -T~ . 10) 

2 V A.7T 2 V ; 
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III. RADIATIVE DECAY OF GRAVITINO 



Radiative decay of gravitino may affect BBN. We analyze this effect assuming that 
gravitinos ip^ mainly decay to photons 7 and photinos 7. 

In order to investigate the photo-dissociation processes, we must know the spectra of the 
high energy photon and electron induced by the gravitino decay. In this section, we will 
derive these spectra by solving the Boltzmann equations numerically. 

Once high energy photons are emitted in the gravitino decay, they induce cascade pro- 
cesses. In order to analyze these processes, we have taken the following radiative processes 
into account. (I) High energy photon with energy e 7 can become e + e~ pair by scattering 
off background photon if the energy of the background photon is larger than m 2 e je 7 ( m e : 
electron mass ). We call this process double photon pair creation. For sufficiently high 
energy photons, this is the dominant process since the cross section or the number density 
of target is much larger than other processes. Numerical calculation shows that this process 
determines the shape of the spectrum of high energy photon for e 7 > m 2 /22T. (II) Below 
the effective threshold of double photon pair creation, high energy photons lose their energy 
by photon-photon scattering. But in the limit of e 7 — > 0, the total cross section for the 
photon-photon scattering is proportional to e 7 and this process loses its significance. Hence 
finally, photons are thermalized by (III) pair creation in nuclei, or (IV) Compton scattering 
off thermal electron. And (V) emitted high energy electrons and positrons lose their energy 
by the inverse Compton scattering off background photon. Furthermore, (VI) the source 
of these cascade processes are the high energy photons emitted in the decay of gravitinos. 
Note that we only consider the decay channel ip^ —> 7 + 7 and hence energy of the incoming 
photon e 7 o is monoclomatic. 

The Boltzmann equations for the photon and electron distribution function f 7 and f e 
are given by 



df-f (e 7 ) df 7 (e 7 ) 



/7V77 

dt 



dt 



+ 



dt Dp 
dt 



+ 



df 7 (e 7 ) 



+ 



cs 



dt pp 

d/ 7 (€ 7 ) 
dt 



+ 



df 7 (e 7 ) 



+ 



ic 



df e (E e ) df e (E e ) 



dt 



+ 



df e (E e ) 



DP 



dt 



+ 



dt 

dt 
df e {E e ) 



PC 



DE 



PC 



dt 



+ 



df e (E e ) 



cs 



dt 



(11) 
(12) 



ic 



where DP (PP, PC, CS, IC, and DE) represents double photon pair creation (photon-photon 
scattering, pair creation in nuclei, Compton scattering, inverse Compton scattering, and the 
contribution from the gravitino decay). Full details are shown in appendix B. 

In order to see the photon spectrum, we have to solve Eq.(ll) and Eq.(12). Since the 
decay rate of gravitino is much smaller than the scattering rates of other processes, gravitinos 
can be regarded as a stationary source of high energy photon at each moment. Therefore, 
we only need a stationary solution of Eq.(ll) and Eq.(12) with non-zero (df 7 /dt)\w at each 
temperature. Note that Eq.(ll) and Eq.(12) are linear equations of f 7 and / e , and hence, 
once Eq.(ll) and Eq.(12) have been solved with some reference value of (d f 7 / <9£)|de we can 
reconstruct the photon spectrum for arbitrary value of (df 7 /dt)\w with T and e 7 o fixed; 
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(df 7 /dt)\ 



/.(e,) = / 7 (c 7 ) x )" J ?'"~"™ . (13) 

For each T and e 7 o, we have calculated the reference spectra / 7 (e 7 ) and f e {E e ) by solving 
Eq.(ll) and Eq.(12) numerically with the condition, 

d/ 7 (e 7 ) = df e (E e ) = 

dt dt { ' 

Typical spectra are shown in Figs.l in which we have shown the case with e 7 o = lOOGeV 
and lOTeV, T = lOOkeV, IkeV, lOeV, and the incoming flux of the high energy photon is 
normalized to be 



d/ 7 (€ 7 ) 



= 6{e 1 - e 7 o) GeV 5 . (15) 

DE 



The behaviors of the photon spectra can be understood in the following way. In the region 
e 7 ^ mg/22T, the photon number density is extremely suppressed since the rate of double 
photon pair creation process is very large. Just below this threshold value, the shape of 
the photon spectrum is determined by the photon-photon scattering process, and if the 
photon energy is sufficiently small, the Compton scattering with the thermal electron is the 
dominant process for photons. In Fig. 2, we have compared our photon spectrum with the 
results of the simple fitting formula used in Ellis et al [16]. As one can see, not only the 
absolute value but also the form of the spectrum differs between them. The fitting formula 
in Ref. [16] is derived from the numerical results given in Refs. [18,19] in which, however, 
the effect of the Compton scattering is not taken into account. Our results indicate that the 
number of Compton scattering events is comparable to that of the inverse Compton events 
for such low energy region, since the number density of the high energy electron is extremely 
smaller than that of high energy photon. Therefore, the deformation of the photon spectrum 
by Compton scattering is expected below the threshold of the photon-photon scattering. 



IV. BBN AND PHOTO-DISSOCIATION OF LIGHT ELEMENTS 

BBN is one of great successes of the standard big bang cosmology. It is believed that 
light elements of mass number less than 7 are produced at early stage of the universe when 
the cosmic temperature is between lMeV and lOkeV. Theoretical predictions for abundances 
of light elements is excellently in good agreement with those expected from observations if 
baryon-to-photon ratio i] B is about 3 x 10~ 10 . 

However the presence of gravitino might destroy this success of BBN. Gravitino may 
have three effects on BBN. First the energy density of gravitino at T ~ lMeV speeds up the 
cosmic expansion and leads to increase the n/p ratio and hence 4 He abundance also increases. 
Second, the radiative decay of gravitino reduces baryon-to-photon ratio and results in too 
baryon-poor universe. Third, the high energy photons emitted in the decay of gravitino 
destroy the light elements. Among three effects, photo-dissociation by high energy photon 
is the most important for gravitino of mass less than ~ ITeV. In the following we consider 
the photo-dissociation of light elements and discuss other effects in Sec. VI. 
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The high energy photons emitted in the decay of gravitinos lose their energy during 
multiple electromagnetic processes described in the previous section. Surviving soft photons 
can destroy the light elements (D, T, 3 He, 4 He) if their energy are greater than the threshold 
of the photo-dissociation reactions. We consider the photo-dissociation reactions listed in 
Table I. For the process D(7,n)p, we have used the cross section in analytic form which is 
given in Ref. [21], and the cross sections for other reactions are taken from the experimental 
data (for references, see Table I). We neglect 4 He(7, D)D and 4 He(7, 2p 2n) since their 
cross sections are small compared with the other reactions. Furthermore, we do not include 
the photo-dissociation process for 7 Li and 7 Be because the cross section data for 7 Be is not 
available and hence we cannot predict the abundance of 7 Li a part of which come from 7 Be. 

The time evolution of the light elements are described by 



n B J2 de i a h^ a ( e i)M e i) + J2 * 7 a^ D (e 7 )n a / 7 (e 7 ), (16) 



= -» He X J E * 7 ^ He ^ a (e 7 )/ 7 (e 7 ) + J2j E rfe 7 a^3 He (e 7 )n a / 7 (e 7 ), (18) 

dm Ee 
dt 



-114 



HeX / * 7 al He ^ a (e 7 )/ 7 (e 7 ) + ^ / de 7 a^4 He (e 7 )n a / 7 (e 7 ), (19) 



where <J l a _^ b is the cross section of the photo-dissociation process i: a+7 — > b+. . . and Ei is the 
threshold energy of reaction %. When the energy of the high energy photon is relatively low, 
i.e. 2MeV ^ e 7 < 20MeV the D, T and 3 He are destroyed and their abundances decrease. 
On the other hand, if the photons have high energy enough to destroy 4 He, it seems that 
such high energy photons only decrease the abundance of all light elements. However since 
D, T and 3 He are produced by the photo-dissociation of 4 He whose abundance is much 
higher than the other elements, their abundances can increase or decrease depending on the 
number density of high energy photon. When the number density of high energy photons 
with energy greater than ~ 20MeV is extremely high, all light elements are destroyed. But 
as the photon density becomes lower, there is some range of the high energy photon density 
at which the overproduction of D, T and 3 He becomes significant. And if the density is 
sufficiently low, high energy photon does not affect the BBN at all. 

From various observations, the primordial abundances of light elements are estimated 
[20] as 

0.22<Y P =( — J < 0.24, (20) 
V Pb J p 

ViA >1.8xl0- 5 , (21) 

p 

< 1.0 x 10~ 4 , (22) 



n H 
n D + n 3H , 



n H 



where /)4 He and pg are the mass densities of 4 He and baryon. The abundances of light 
elements modified by gravitino decay must satisfy the observational constraints above. In 
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order to make precise predictions for the abundances of light elements, the evolutional 
equations (16) - (19) should be incorporated with the nuclear network calculation of BBN. 
Therefore, we have modified Kawano's computer code [28] to include the photo-dissociation 
processes. 

From the above arguments it is clear that there are at least three free parameters, i.e. 
mass of gravitino 7713/2, reheating temperature Tr and baryon-to-photon ratio i] B . Further- 
more we also study the case in which gravitino has other decay channels. In the present 
paper we do not specify other decay channel. Instead, we introduce another free parameter 
B 1 which is the branching ratio for the channel ip^ — > 7 + 7. Therefore we must study the 
effect of gravitino decay on BBN in four dimensional parameter space. However in the next 
section it will be shown that the baryon-to-photon ratio r] B is not important parameter in 
the present calculation because the allowed value for r] B is almost the same as that in the 
standard case (i.e. without gravitino). 

V. RESULTS 

A. B 1 — 1 case 

First we have investigated the photo-dissociation effect when all gravitinos decay into 
photons and photinos {B 1 — 1). We have taken the range of three free parameters as 
lOGeV < m 3 / 2 < lOTeV, 10 5 GeV < T R < 10 13 GeV and 1(T 10 < f] B < 10~ 9 . In this 
calculation, we assume that photino is massless. The contours for the critical abundances 
of the light elements D, (D+ 3 He) and 4 He in the r] B — Tr plane are shown in Figs. 3 for 
(a)m 3 / 2 = lOGeV, (b)100GeV, (c)lTeV and (d)10TeV, respectively. For low reheating 
temperature {Tr < 10 6 GeV), the number density of the gravitino is very low and hence the 
number density of the induced high energy photons is too low to affect the BBN. Therefore 
the resultant abundances of light elements are the same as those in the standard BBN. The 
effect of the photo-dissociation due to gravitino decay becomes significant as the reheating 
temperature increases. 

As seen in Figs. 3, the allowed range of baryon-to-photon ratio is almost same as that 
without gravitino for 777,3/2 ^5 ITeV, i.e. very narrow range around r] B ~ 3 x 10~ 10 is allowed. 
However for 7773/2 ~ ITeV and Tr ~ 10 9 GeV or 7773/2 ~ ITeV and Tr ~ 10 12 GeV , lower 
values of r] B are allowed (Fig. 3(c)). In this case, the critical photon energy (~ 777 2 /22T) 
for double photon pair creation process is lower than the threshold of photo-dissociation 
reaction of 4 He. Therefore, for Tr < 10 12 GeV, the abundance of 4 He is not affected by the 
gravitino decay. Then the emitted photons only destroy 3 He and D whose abundances would 
be larger than the observational constraints for low baryon density if gravitino did not exist. 
Therefore one sees the narrow allowed band at T — 10 9 GeV where only a small number of 
3 He and D are destroyed to satisfy the constraints (21) and (22). For Tr 10 12 GeV, since 
a large number of high energy photons are produced even above the threshold of double 
photon pair creation, a part of 4 He are destroyed to produce 3 He and D, which leads to 
the very narrow allowed region at Tr ~ 10 12 GeV. However even in this special case, the 
upper limit of allowed reheating temperature is changed very little between r] B — 10~ 10 and 
r\B ~ 3 x 10~ 10 . This allows us to fix T] B — 3.0 x 10~ 10 in deriving the upperbound of the 
reheating temperature. 
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The allowed regions that satisfy the observational constraints (20)-(22) also shown in 
Figs. 4 in the — Tr plane for rjg — 3 X 10~ 10 . In Figs. 3 and Fig. 4(a) one can see four 

typical cases depending on T# and m^/2- 

• m 3/2 < ITeV, T R < 10 n GeV: 

In this case the lifetime of the gravitino is so long that the critical energy for double 
photon process (~ m^/22T) is higher than the threshold of the photo-dissociation 
reactions for 4 He at the decay time of gravitino. Thus 4 He is destroyed to produce T, 
3 He and D. (Since T becomes 3 He by /?-decay, hereafter we mean T and 3 He by the 
word " 3 He".) Since the reheating temperature is not so high, the number density of 
gravitino is not high enough to destroy all the light elements completely. As a result, 
3 He and D are produced too much and the abundance of 4 He decreases. To avoid 
the overproduction of 3 He and D, the reheating temperature should be less than ~ 
(10 6 - 10 9 )GeV. 

• m 3/2 < ITeV, T R > 10 n GeV: 

The lifetime is long enough to destroy 4 He and the gravitino abundance is so large 
that all the light elements are destroyed since the reheating temperature is extremely 
high. This parameter region is strongly excluded by the observation. 

• ITeV < m 3 / 2 £ 3TeV: 

The lifetime becomes shorter as the mass of gravitino increases, and the decay occurs 
when double photon pair creation process works well. If the cosmic temperature at 
t — T3/2 is greater than ~ mg/22i?4 He (where i?4 He ~ 20MeV represents the typical 
threshold energy of of 4 He destruction processes), 4 He abundance is almost unaffected 
by the high energy photons as can be seen in Fig. 3(c). In this parameter region, 
overproduction of (D+ 3 He) cannot occur since 4 He is not destroyed. In this case, the 
destruction of D is the most important to set the limit of the reheating temperature. 
This gives the constraint of T R < 10 9 - 10 12 GeV. 

3TeV: 

In this case the decay occurs so early that all high energy photons are quickly ther- 
malized by double photon process before they destroy the light elements. Therefore 
the effect on BBN is negligible. Fig. 3(d) is an example of this case. The resultant 
contours for abundances of light elements are almost identical as those without the 
decay of gravitino. 



B. B 1 < 1 case 

So far we have assumed that all gravitinos decay into photons and photinos. But if other 
superpartners are lighter than gravitino, the decay channels of gravitino increases and the 
branching ratio for the channel ip^ — > 7 + 7 becomes less than 1. In this case, various decay 
products affect the evolution of the universe and BBN. In this paper, instead of studying all 
decay channels, we consider only the 7 + 7 channel with taking the branching ratio B 1 as an- 
other free parameter. With this simplification, the effect of all possible decay products other 
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than photon is not taken into account. Therefore the resultant constraints on the reheat- 
ing temperature and the mass of gravitino should be taken as the conservative constraints 
since other decay products may destroy more light elements and make the constraints more 
stringent. 

Although we have four free parameters in the present case, the result for B 1 — 1 implies 
that the allowed range of Tr and 7773/2 is obtained if we take baryon-to-photon ratio to be 
•3 x 10~ 10 . Since our main concern is to set the constraints on Tr and 7713/2, we can safely 
fix T] B (= 3 x 10" 10 ). 

The constraints for B 1 — 0.1 and B 1 — 0.01 is shown in Fig. 4(b) and Fig. 4(c) which 
should be compared with Fig. 4(a) (B 1 — 1 case). Since the number density of the high 
energy photons is proportional to J3 7 , the constraint on the reheating temperature becomes 
less stringent as B 1 decreases. In addition, the total lifetime of gravitino is given by 

t 3 /2 = t(W^7 + 7)x5 7 . (23) 

Thus the gravitinos decay earlier than that for B 1 — 1 case and the constraints from ( 3 He 
+ D) overproduction becomes less stringent. This effect can be seen in Fig. 4(b), where the 
constraint due to the overproduction of ( 3 He + D) has a cut at 777,3/2 — 400GeV compared 
with ~ ITeV for B 1 = 1. 

In Fig. 5, the contours for the upperbound of reheating temperature are shown in the 
"73/2 — B 1 plane. One can see that the stringent constraint on Tr is imposed for 777,3/2 ^5 
lOOGeV even if the branching ratio is small. As mentioned before this constraint should be 
regarded as the conservative one and the actual constraint may become more stringent by 
the effect of other decay products, which will be investigated elsewhere. 



VI. OTHER CONSTRAINTS 

In the previous section, we have considered the constraints from the photo-dissociation 
of light elements. But as we have seen, if the mass of gravitino is larger than a few TeV, 
gravitino decay does not induce light element photo-dissociation and no constraints has 
been obtained. In the case of such a large gravitino mass, we must consider other effects of 
gravitino. 

If we consider the present mass density of photinos produced by the gravitino decay, we 
can get the upperbound of the reheating temperature. In SUSY models with i?-invariance 
(which is the usual assumption), the lightest superparticle (in our case, photino) is sta- 
ble. Thus the photinos produced by the decay of gravitinos survive until today, and they 
contribute to the energy density of the present universe. Since one gravitino produces one 
photino, we can get the present number density of the photino: 

n 7 = y 3/2 (T < lMeV) x ^T 3 , (24) 

where T$ is the present temperature of the universe. Density parameter of photino il 7 = 
777 7 77 7 //7 c can be easily calculated, where is the photino mass, /9 C ~ 8.1 x 10 _47 /i 2 GeV 4 is 
the critical density of the universe and h is the Hubble parameter in units of lOOkm/sec/Mpc. 
If we constrain that < 1 in order not to overdose the universe, the upperbound of the 
reheating temperature is given by 
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T ^ 2 ' 7xioii (Toi^r ,i2Gev - < 25 > 

where we have ignored the logarithmic correction term of T, tot . To set the upperbound 
on the reheating temperature, we need to know the mass of photino. If one assumes the 
gaugino-mass unification condition, the lower limit on the mass of photino is 18.4GeV [29]. 
Then we can get the following upperbound on the reheating temperature: 

T R < 7.1 x 10 n /i 2 GeV. (26) 

Note that this bound is independent of the gravitino mass and branching ratio. 

Another important constraint comes from the effect on the cosmic expansion at BBN. 
As mentioned before, if the density of gravitino at nucleosynthesis epochs becomes high, the 
expansion of the universe increases, which leads to more abundance of 4 He. We study this 
effect by using modified Kawano's code and show the result in Fig. 6. In the calculation, we 
take t]b — 2.8 x 10~ 10 and r n — 887 sec (where r n — (889±2.1)sec is the neutron lifetime [34]) 
so that the the predicted 4 He abundance is minimized without conflicting the observational 
constraints for other light elements. The resultant upperbound of the reheating temperature 
is given by 

^<2xl0 13 Gev(^)~\ (27) 

for m3/2 > ITeV. 

VII. DISCUSSIONS 

We have investigated the production and the decay of gravitino, in particular, the effect 
on BBN by high energy photons produced in the decay. We have found that the stringent 
constraints on reheating temperature and mass of gravitino. 

Let us compare our result with those in other literatures. Our number density of gravitino 
produced in the reheating epochs of the inflationary universe is about four times larger than 
that in Ellis et al [11]. Since Ellis et al [11] note that they have neglected the interaction 
terms between gravitino and chiral multiplets (which is the second term in Eq.(Al)), they 
might underestimate the total cross section for the production of gravitino. All previous 
works concerning gravitino problem were based on the gravitino number density given by 
Ellis et al [11]. Therefore our constraints are more stringent than others. Furthermore, our 
photon spectrum is different from that in Ref. [16] as shown in Fig. 2. The spectrum adopted 
by Ref. [16] has more power to destroy light elements above threshold for the photon-photon 
scattering and less power below the threshold. It is expected that the difference comes mainly 
from the neglect of Compton scattering in Ref. [16]. In addition, the absolute value of our 
spectrum slightly depends on the energy of initial photons. This can be explained in the 
following way. If the injected photons become more energetic, they induces larger number 
of scattering events before thermalized, which produce larger number of secondary high 
energy photons. Therefore, the energy of the initial photon gets larger, the number density 
of high energy photons becomes higher. This effect was not taken into account in Ref. [16]. 
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Therefore we believe that the spectrum that we have obtained is more precise than those 
used in other works. 

In summary, we have investigated the photo-dissociation processes of light elements due 
to the high energy photons emitted in the decay of gravitino and set the upperbound of the 
reheating temperature by using precise production rate of gravitino and the spectrum of 
high energy photon. Together with other constraints ( the present mass density of photino 
and enhancement of cosmic expansion due to gravitino) we have obtained the following 
constraint; 

m 3/2 < lOOGeV, (28) 
lOOGeV < m 3/2 < ITeV, (29) 
ITeV < m 3/2 < 3TeV, (30) 
3TeV < m 3/2 < lOTeV. (31) 

This provides a severe constraint in building the inflation models based on supergravity. 
In this paper we also study the gravitino which decays into other channels by taking the 
branching ratio as a free parameter. Although this gives conservative upperbound for the 
reheating temperature, the precise constraints cannot obtained unless various processes in- 
duced by other decay products are fully taken into account. This will be done in the future 
work [35]. 



T R < 10 6_7 GeV 
T R < 10 7 " 9 GeV 
T R < 10 9_12 GeV 
T R < 10 12 GeV 
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APPENDIX A: INTERACTION OF GRAVITINO 

In this appendix, we will discuss the interaction of the gravitino with the ordinary par- 
ticles. From the supergravity lagrangian [5], we can obtain the relevant interaction terms of 
gravitino ip^ with gauge multiplets (A fl7 A) and chiral multiplets (0, %); 

£ = ,s.'l/ /V/ ' ^ 7 " 7 ^ ^ Fvp + { \/2.l/' V/ " "" /,V/ /; "" + k ' C '} ' ^ 

Note that other interaction terms including gravitino field are not important for our analysis 
since their contributions are suppressed by factor M _1 . 

Combining Eq.(Al) with the renormalizable part of the SUSY lagrangian, we have cal- 
culated the helicity ±| gravitino production cross sections and the results are shown in 
Table II. Note that the cross sections for the processes (B), (F), (G) and (H) are singular 
because of the t-channel exchange of gauge bosons. These singularities should be cut off 
when the effective gauge boson mass m e ff due to the plasma effect has been taken into 
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account. Following Ref. [11], we take 6 = (1 =F cos #)mm = ( m ejf/2^ 2 ) where 9 is a scat- 
tering angle in the center of mass frame, and in our numerical calculations, we have chosen 
log(m 2 ejJ /T 2 ) = 0. 

The total cross section in thermal bath T, tot is defined by 

^tot — 7> XI ^^Vy ^(x+y^^+z) , (A2) 
x,y,z 

where <j( x+y ^ +z ) is the cross section for the process x + y —>■ ip^ + z, i] x — 1 for incoming 
bosons, r] x — \ for fermions. For the MSSM particle content, E tot is given by 

^ot = {2-50^(T) + 4.9% 2 2 (T) + 11.78<7 3 2 (T)} , (A3) 

where c/i, 52 an d 93 are the gauge coupling constants of the gauge group U(l) y , SU(2) L 
and SU(3) C , respectively. Note that in high energy scattering processes, effect of the renor- 
malization group flow of the gauge coupling constants should be considered. Using the one 
loop /^-function of MSSM, solution of the renormalization group equation of gauge coupling 
constants is given by 

^^-^-^MwX' 2 ■ (A4) 

with bi — 11, 62 = 1, &3 = — 3. In this paper, we have used the above formula. 

From Eq.(Al), we can also get the decay rate of the gravitino. In this paper, we only 
consider ip^ — > 7 + 7 decay mode, for which the decay rate is given by 

where is the photino mass. In the case of <C ^3/2, this decay rate corresponds to 
the lifetime 

t{% - 7 + 7) = 3.9 x 10 8 (j^g^) 3 sec. (A6) 




APPENDIX B: BOLTZMANN EQUATION 

In order to calculate the high energy photon spectrum, we must estimate the cascade 
processes induced by the radiative decay of the gravitinos. In our calculation, we have 
taken the following processes into account; (I) double photon pair creation, (II) photon- 
photon scattering, (III) pair creation in nuclei, (IV) Compton scattering off thermal electron, 
(V) inverse Compton scattering off background photon, and (VI) radiative decay of the 
gravitinos. The Boltzmann equations for this cascade processes are given by 
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Below, we see contributions from each processes in detail. 



(I) DOUBLE PHOTON PAIR CREATION [ 7 + 7 -»■ e+ + e" 7 

For the high energy photon whose energy is larger than ~ m^/22T, double photon pair 
creation is the most dominant process. 

The total cross section for the double photon pair creation process a dp is given by 



a DP (P) = \vrl (l - f) I (3 - ±f) log \ 



P 



(B3) 



where r e — a/m e is the classical radius of electron and (3 is the electron ( or positron ) 
velocity in the center of mass frame. Using this formula, one can write down (a/ 7 /at)|op as 

2m! 



dfJe~) 



dt 



/ 7 ( e 7 ) 



ic 



c 2 



'i -2 

'e/ c 7 



ds (T DP (f3) 



(3=^/l-(l/s) 



(B4) 



The spectrum of the final state electron and positron is obtained in Ref. [30], and 
(a/ e /at) I dp is given by 



df e (E e ) 
dt 



x 2 4 
= -7rr„m„ 



DP 



de 



7 —2 V 77 t 7/? 



(B5) 



where / 7 represents the distribution function of the background photon at temperature T, 



x 



1 



7r 2 exp(e 7 /T) — 1 



(B6) 



and function G(E e , e 7 , e 7 ) is given by 



v\ 2 



G^,^,^) = v : ^ e; log :?r e \„ -8-^ 



E P E' 



mUE e + E' e ) 



mi 



2{2e 7 (E e + E' e )-mi\ (E e + E[ 



m 2 e E e E' e 



E e + E f 



e 7 (E e + E' e ) E*E' e 2 



(B7) 



with 



E p — e 7 + e 7 



E P 
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(II) PHOTON-PHOTON SCATTERING ^ + 7^7 + 7 ] 



If the photon energy is below the effective threshold of the double photon pair creation, 
photon-photon scattering process becomes significant. This process is analyzed in Ref. [19] 
and for e' 7 < 0(m 2 e /T), (df~ t /dt)\ PP is given by 



dm { ' 1 



dt 



2 



35584 „J!„4 f(r V2 J 1 S , / e 7 



pp 10125^™" k dc ^W \ 1 " % + [t ) > /• 



oc 



1946 
10125tt 



/7(s)« r e m e 6 S 3 / de 7 e^/ 7 (e 7 ). (B8) 



For a larger value of e we cannot use this formula. But in this energy region, photon-photon 
scattering is not significant because double photon pair creation determines the shape of the 
photon spectrum. Therefore, instead of using the exact formula, we have taken m 2 /T as a 
cutoff scale of (df 7 /dt)\ PP7 i.e., for e 7 < fn 2 e /T we have used Eq.(B8) and for e 7 > m 2 /T we 
have taken 



a/ 7 ( C ; > m 2 jT) 



= 0. (B9) 

pp 



dt 

Note that we have checked the cutoff dependence of spectra is negligible. 



(Ill) PAIR CREATION IN NUCLEI [ 7 + N -»■ e+ + e" J 



Scattering off the electric field around nucleon, high energy photon can produce electron 
positron pair if the photon energy is larger than 2m e . Denoting total cross section of this 
process o P c 1 (9/ 7 /9£)|np is given by 



dt 



-n N apcf 7 (e 7 ) 7 



(BIO) 



NP 



where is the nucleon number density. For <Jpc 7 we have used the approximate formula 
derived by Maximon [31]. 

Differential cross section for this process dupc/dE e is given in Ref. [32], and {df e /dt)\^ P 
is given by 



df e (E e ) 



dt 



= n N 



NP 



E c +m t 



, dope , , x 
7 Ji\ e i)' 



(Bll) 



f/Fj COMPTON SCATTERING [ 7 + e" -»■ 7 + e" 7 

Compton scattering is one of the processes by which high energy photons lose their 
energy. Since the photo-dissociation of light elements occurs when the temperature drops 
below ~ O.lMeV, we can consider the thermal electrons to be almost at rest. Using the total 
and differential cross section at the electron rest frame <Jcs an d docs I dE e ^ one can derive 
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where n e is the number density of thermal electron. 



(V) INVERSE COMPTON SCATTERING [e ± +j^e ± +j ] 



Contribution from the inverse Compton process is given by Jones [33], and (df/dt)\ic is 
given by 
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dt 

dfe(E' e , 
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= Znr e m e 
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where function F(e 7 , i? e , e 7 ) is given by 



r £ ?) s 



F(e 7 ,£ e , e 7 ) = 2glogg + (1 + 2g)(l - q) + - ± ™ ^ (1 - g) 



(B16) 



with 



r _ 4e 7 E e 



(T/j GRAVITINO RADIATIVE DECAY / W -»■ 7 + 7 7 

Source of the non-thermal photon and electron spectra is radiative decay of gravitino. 
Since gravitinos are almost at rest when they decay and we only consider two body decay 
process, incoming high energy photons have unique energy e 7 o, which is given by 



e 7 o — 



m i/2 - ™% 



2m 



3/2 



(B17) 



Therefore, (9/ 7 /9£)|de can be written as 
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FIGURES 
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FIG. 1. Typical spectra of photon (the solid lines) and electron (the dotted lines). We have 
taken the temperature of the background photon to be T — lOOkeV, IkeV, lOeV, and the energy 
of the incoming high energy photon e 7 o is (a) lOOGeV and (b) lOTeV. Normalization of the initial 
photon is given by e 7 o x (9/ 7 (c 7 )/9<)|de — ^( e 7 — e 7o) GeV 5 . 
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FIG. 2. Photon spectrum derived from the fitting formula used in Ref. [16] is compared with 
our results. We have taken the temperature of the background photon to be lOOeV and the 
normalization of the incoming flux is the same as Fig.l. The solid line is the result of fitting 
formula, and the rest of two are our results with e 7 = lOOGeV (dotted line) and e 7 = lOTeV 
(dashed line). 
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FIG. 3. Contours for critical abundance of light elements in the rjg — Tr plane for (a) 
m 3/2 = lOGeV, (b) m 3/2 = lOOGeV, (c) m 3/2 = ITeV and (d)m 3/2 = lOTeV. 
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FIG. 4. Allowed regions in m^j2 — Tr plane for (a) B 1 — 1, (b)i? 7 = 0.1 and (c) B 1 — 0.01. In 
the region above the solid curve 3 He and D are overproduced, the abundance of 4 He is less than 
0.22 above the dotted curve and the abundance of D is less than 1.8 x 10 above the dashed curve. 
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FIG. 5. Contours for the upper limits of the reheating temperature in the 771.3/2 — -B7 plane. 
The numbers in the figure denote the limit of the reheating temperature. 
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FIG. 6. Upperbound of the reheating temperature. Dashed line represents the constraint from 
the present mass density of photino. Solid curve represents the upperbound requiring 4 He < 0.24. 
Constraints from D photo-dissociation is also shown by dotted line. 
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TABLES 



Reaction 


Threshold (MeV) 


References 


D +7 — y n + p 


2.225 


[21] 


T +7 — ► n+ D 


6.257 


[22], [23] 


T +7 — y p + n + n 


8.482 


[23] 


3 He +7 -yp+D 


5.494 


[24] 


3 He +7 -yp+D 


7.718 


[24] 


4 He +j^p+T 


19.815 


[25] 


4 He +7 -»• n+ 3 He 


20.578 


[26] 


4 He +i —yp + n+ D 


26.072 


[27] 



TABLE I. Photo-disociation reactions 



Process 

(A) 

(B) 

(C) 

(D) 

(E) 

(F) 
(G) 

(H) 
(I) 
(J) 



A a + A b 

A a + A 6 • 

A a + <f>i - 

Xi + <l>j - 
X a + A 6 - 

A a + Xt - 
A a + <f>i - 

Xi + Xj - 

4>i + <t>) - 



y-ip + \ c 
>ip + A c 

>^ + Xj 
ytp + (f)j 
tp + A a 
► ip + X c 

■■<P + Xj 
■tp + (f>j 
■ip + X a 
ip + \ a 



a = (c, 2 /64vrM 2 )x 



(8/3) f abc 
h 2 

^ jaoc 



rpa 

ji 

rpa 

ji 

rpa 

ji 



{-(3/2) + 21og(2/6) + 6-(l/8)6 2 ] 



jabc {-(62/3) + 161og[(2 - 8)/ 8] + 228 - 28 2 + (2/3)8 3 } 



ji 
2 



3 l 

(8/3) 
(16/3) 



{-2 + 21og(2/6) + 6} 
{-6 + 81og(2/6) + 46- (l/2)6 2 } 

2 

2 



J* 



rpa 

ji 



TABLE II. Total cross sections for the helicity ±| gravitino production process. Spins of the 



initial states are averaged and those of the final states are summed, f and T-j represent the 
structure constants and the generator of the gauge group, respectively. Note that for the processes 
(B), (F), (G) and (H), we cut off the singularities due to the t-, li-channel exchange of gauge bosons, 
taking (1 ± cos 0^ m i n — 6 where is the scattering angle m the center-of-niass frame. 
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